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CHAPTER I 


INTRODUCTION 

In the design and use of large composite panels in structures such 
as aircraft and space vehicles a major concern is the ability of the 
panel to function (continue to carry a substantial part of its design 
load) after being damaged. One technique that has evolved from the 
stringer reinforced metallic panel is the buffer strip or hybrid panel 
shown in Figure 1. Because of the fabrication methods used in composites 
it is possible to make such a laminate by replacing specific fibers, 
usually parallel to the load axes, with fibers of the appropriate physical 
and geometric properties necessary to arrest a crack that originates in 
the parent laminate material. Since the buffer strips are usually narrow 
and relatively far apart the stiffness, weight, and strength of the 
undamaged laminate is not significantly effected by the replacement. 

Much experimental work has been done to investigate this behavior, 
i.e. to determine the best buffer strip material, with the studies of 
Eisenmann and Kaminski [1], Hess, Huang and Rubin ,[2J , Avery and porter 13] 
Verette and Labor [4], and Poe and Kennedy [5] being significant contribu- 
tions in this area. The same cannot be said for published analytical 
solutions and it is this question that is considered in the present study. 

Some of the first work in modeling a uni -directional composite con- 
taining broken fibers was presented by Hedgepeth [6] where the case of 
no additional damage other than the initial notch was considered. This 
study was extended by Hedgepeth and Van Dyke for the special case of one 


main laminate 




Figure 1. A Typical Buffer Strip Laminate 





broken fiber with matrix yielding parallel to the fiber [7] and for one 
fiber with longitudinal splitting in the matrix [8]. Goree and Gross 

[9] extended the Hedgepeth solutions to include longitudinal matrix 
yielding and splitting for an arbitrary number of broken fibers. The 
results of Goree and Gross gave very good agreement with experimental 
results for brittle matrix composites which exhibit large longitudinal 
matrix splitting. For ductile matrix composites such as boron/aluminum, 
which exhibit large yielding but very little splitting in the matrix, 
this model predicted the right trend but the agreement was not very 
good, especially for short notch lengths. Goree, Dharani and Jones 

[10] attempted two modifications to the above solution. First, the 
matrix was assumed to be strain-hardening and secondly, a cover sheet 
was included over the main laminate. The results of [10] showed that 
the inclusion of either a strain-hardening matrix or the addition of a 
cover sheet did not improve the agreement between the predicted and 
experimental results. Based on the observed fact that in addition to 
longitudinal yielding of the matrix, a certain amount of stable trans- 
verse extension of the initial notch under increasing applied load 
takes place, Dharani, Jones and Goree [11] then extended the solution of 
[9] to include transverse damage ahead of the initial notch in addition 
to the longitudinal matrix damage. The results of [11] showed a very 
significant improvement in the ability of the model to represent the 
behavior of a ductile matrix composite. 

In all these analytical studies the laminate is modeled as a two- 
dimensional region having a single row (mono-layer) of parallel, identi- 
cal, equally spaced fibers, separated by matrix. The damage is taken to 
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consist of an arbitrary number of broken fibers such that all breaks 
lie along the x-axis, but they need not form a continuous break (notch). 
The fibers are assumed to be of much higher strength and extensional 
stiffness than the matrix and all the axial load is assumed to be 
carried by the fibers, with the matrix transfering load by shear 
stresses as given by the classical shear-lag assumption. One very impor- 
tant feature of the shear-lag assumption is that is simplifies the equili- 
brium equations by removing the transverse displacement dependence from 
the longitudinal equilibrium equation. The fiber stress and matrix shear 
stress can then be determined without solving the transverse equation. 

The methods of analysis developed and discussed in the above 
studies [6-11] are extended in the present work to determine the 
stresses and displacements in a hybrid uni -directional laminate having 
an initial notch in the vicinity of a single finite width region of 
different material properties as shown in Figure 2. This geometry is 
an idealization of the usual periodic placement of buffer strips as 
indicated in Figure 1 and it is assumed that the stresses near the 
notch tip and the single buffer strip are approximately the same as 
those in a wide panel with relatively narrow buffer strips. 

Of particular interest is the investigation of the behavior of the 
laminate as a function of the relative ultimate stress and extensional 
modulus of the buffer strip fibers, the buffer strip width and initial 
notch length. 

As an initial step in understanding the basic mechanism of crack 
growth and arrest in a hybrid laminate, the main thrust of the investiga- 
tion will be to study the behavior of the laminate as a function of the 
design parameters; fiber/matrix properties, buffer strip geometry and 
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the initial crack length. In order to simplify the analysis, the 
effect of any additional damage, e.g. longitudinal matrix yielding 
and splitting [9] or transverse matrix and fiber damage [11] will not 
be considered. 

A typical buffer strip laminate usually contains angle plies as 
well as zero degree plies. It is felt, however, that much of the 
characteristic behavior of the buffer strip region can be represented 
by the uni-directional laminate, as a major portion of the load is 
carried by these fibers. It appears that a primary function of the 
angle plies is to prevent longitudinal matrix splitting in a brittle 
matrix such as epoxy. This can be accounted for to some degree in the 
present solution by allowing the matrix to support large strains without 
spl i tting. 

The presentation of the solution will follow the order of the 
development by the author as this seems to indicate more clearly the 
significant points of the analysis. Results will, however, only be given 
for the final solutions corresponding to the buffer strip laminate and 
the finite width strip. First the solution for a crack in a half-plane 
having arbitrary shear stresses applied to the free-edge, (which forms 
the basis for all the later solutions), will be developed. By matching 
boundary conditions along the interface, the solution for two different 
adjoined half-planes will be given and then the adjoined half-plane solu- 
tion will be modified to account for a second interface, resulting in the 
buffer strip problem. The case of the finite width strip with broken 
fibers is obtained as a special case of the buffer strip problem by set- 
ting the shear stresses along the two interfaces to zero. These solutions 
are all presented in Chapter II. Chapter III deals with the numerical 
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technique used to solve a system of linear algebraic simultaneous equa- 
tions coupled with a set of linear Fredholm integral equations of the 


second kind. A special method developed to evaluate accurately inte- 
grals having a cusp is emphasized. Chapter IV gives results and con- 
clusions for the buffer strip problem and the finite width strip. 
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CHAPTER II 


FORMULATION 

Uni-Direct.-jonal Half-Plane with Broken Fibers 

A uni -directional array of parallel fibers with an arbitrary number 
of broken fibers in the form of a notch is shown in Figure 3. The 
laminate is subjected to a prescribed shear stress, x,(y), along the 

a 

free edge in addition to a remote uniform tensile stress in the axial 
direction. Fiber breaks occur along the x-axis (axis of symmetry) and, 
since the loading is symmetric, only the upper half of the laminate is 
considered in the analysis. 

The fibers are taken to be of much higher strength and extensional 
stiffness than the matrix and therefore all of the axial load is assumed 
to be carried by the fibers with the matrix transferring load by shear 
stresses as given by the classical shear-lag assumption. The axial fiber 
stress, i matrix shear stress, x^(y), are then given by the 

simple relations 

civ (y) 

t„(y) = -jf [v„(y) - v^_.,{y)] . (1) 

Where v^(y) is the axial displacement of the fiber n at the loca- 
tion y, Ep is the Young's modulus of the fiber, Gj^j is the equivalent 
matrix shear modulus and h is a shear transfer distance. Because of the 
interference between fibers it is unlikely that 6j^ will be the homo- 
geneous matrix shear modulus or h the actual fiber spacing, and it is 
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Figure 3. Uni-Directional Half-Plane with Broken Fibers 





pointed out in [10] that these values should be determined experimentally 
for a given laminate. It is also shown in [10] that a single experi- 
ment giving the crack opening displacement as a function of applied 
load is sufficient to determine both the equivalent shear modulus G|y| 
and the shear transfer distance h for a particular laminate, indepen- 
dent of the notch length. That is, these parameters are material con- 
stants and depend only on the fiber and matrix properties and the 
fiber volume fraction but not on the sizes of the damage. 

By virtue of the shear-lag assumption the longitudinal and trans- 
verse equilibrium equations become decoupled and the fiber axial dis- 
placements and stresses can be obtained without solving the transverse 
equilibrium equation. Therefore, only the equilibrium equation in the 
longitudinal (axial) direction will be considered. With reference to 
the free-body diagram of a typical fiber-matrix region shown in Figure 3, 
the equilibrium equations in the longitudinal direction is given by 

Ap dcr (y) 

1 ° ’ ( 2 ) 

for all fibers except n = 0, and 

Ap da^(y) 

"t" — dy — ^ ^ fiber 0 . (3) 

Using the stress-displacement relations. Equation (1), in the above 
equilibrium equations, the following set of differential -difference equa- 
tions is obtained: 

ApEph d^v 

* Vi - + v^.i = 0 , and (4) 

ApEph d^v 

_ = . (5) 
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Noting the coefficient of the second derivative term in the above equa- 
tions, the following changes in the variables are suggested: 


let y = 


1 


ApEph 




n , 


dv. 


» — P* ■ • 

(7 = 0* a = tr- “ji — « 

n “ n F dy ’ 


and 


Aph 


Algebraic manipulation then gives 


( 6 ) 


dV 

= a —j — and = a 
n “ dn n ° 


1 




M^F 


Epht 


(Vn-Vl) • 


(7) 


where, n » and V^(n) are non-dimensional. 


The resulting equilibrium equations in non-dimensional form are given by 


d^V. 


dn 


+ V - 2V + V T = 0 , and 
2 n+1 n n-1 ’ 


d^V 


dn' 




(8) 


(9) 


where. 


yn) = 


Epth T^(y) 

Ar-G|y, a 
F M “ 


These differential -difference equations are reduced to differential 
equations by introducing the even valued transform as 


V(n,e) = I V (n) cos[(n + i)e] , 
n=0 

from which 

Vj^(n) = f / V(n,e)cos[(n + ^)e]d0 . 


(10) 


(11) 


Making use of the above transformation and the orthogonality property of 
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the circular functions, the two equilibrium equations may be written as 
one equation valid for all values of n and n as follows: 

— / 2[1- cos( 0 )]vl cos[(n + i) 0 ]de 

"" 0 Idn^ J ^ 

= - / T^(n)cos(0/2)cos[(n + i) 0 ]d 0 , 

IT 0 a 

which is of the form 


o IT 1 

— / F(rij9)cos[(n + -«)0]d0 = 0 for all n and n . 

^0 

Noting the definition of V(ri 50 ) in Equation (10) and (11) it is seen that 
the function F(n»0) is even valued in 0 and therefore, if the integral is 
to vanish for all ns the function F(ns0) must be zero. The single equa- 
tion specifying V(ns0) is then 


2 - 

- 6^V = T^(n)cos(0/2) , 
dn 

where, 6^ = 2[1- cos(0)] = 4 sin^(0/2) . 


( 12 ) 


The solution to the problem of vanishing stresses and displacements 
at infinity and uniform compression on the ends of the broken fibers will 
now be sought. The complete solution is obtained by adding the results 
corresponding to uniform axial stress and no broken fibers to the follow- 
ing solution. The appropriate boundary conditions are as follows: 
dV_(n) 

Vj^(n) = 0 , — — = 0 , as n ^ for all fibers, (13) 

- 

— ^ — = a^(n) = ^- 1 , at n = 0 , for all broken fibers, (14) 

V^(n) = 0 , at n = 0 , for all unbroken fibers. (15) 

The complete solution to Equation (12), satisfying vanishing stresses 
and displacements at infinity, is given by 
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( 16 ) 


V(n,e) = A(e)e"®'^ - / s1nh[6(n-t)]^g(t)dt , 

n 

where the function A(e) is yet unknown. The remaining two boundary con- 
ditions give 

dV (0) 2 ^ ” - 1 

— — = — / [-<5A(e) - cos(e/2)/ cosh(6t)T^(t)dt]cos[(n + -^)e]d0 = -1 , 

" ° ° (17) 

for all broken fibers, and 


V (0) = — / [A(e) + f sinh(6t)T_(t)dt]cos[(n + i)e]d 0 = 0 , 


(18) 


for all unbroken fibers. Equation (18) is solved exactly by taking 


A(0) + 


M 


W I’J ■! 

f sinh(6t)x_(t)dt = 2 B_ cos[(N* + m + i-)0] , (19) 

m c 

0 m=l 


where, M is the number of broken fibers. By eliminating A(e) between 
equations (17) and (19), the stress boundary condition reduces to 


o IT M 1 .. 

- / £ B^ cos[(N* + m + -j)0]cos[(n + :^)0]6d0 

0 m=l ^ 

2 7T 1 “ 

+ — / cos(e/2)cos[(n + -^)e] / e x^(t)dt d© = 1 , (20) 

o 0 

for n = N* + 1 , . . . ,N. 

For a given shear stress distribution, x_(t). Equation (20) reduces 

a 

to a set of linear algebraic equations in the Fourier constants B^^ . From 
Equations (16) and (19), A(e) may be eliminated to obtain V(n»0) in terms 
bf the constants B^ . Recalling the relation between V(ri,0) and V^(n)» 
an expression can be obtained for the axial fiber displacement V^(n) as 

o 7T M 11 

V„(n) = — / e 2 ! B cos[(N* + m + i)0]cos[(n + i)o]d0 

" 0 m=l ^ 

- - / / D(6,n,t) ^^(t)dt cos[(n+l)0]d0 , (21) 

o 0 

where D(6,ri,t) = e - e . 
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The axial fiber stress is obtained by differentiating Equation (21) with 

respect to n and is given by 

dV^(n) 

^n dn 

o TT -(5n M 1 -t 

= - ^ / 6e s cos[(N* + m + ^)e]cos[(n + Y) 0 ]de 

o m=l 


■1 77 CO f -6(TV^"t) ~6 ri"t 

I/. __ 


■1 7T «> f 

- — / cos(0/2)/ { e 
^0 0 


- P e 


where, p = 1 for t £ n » and 


= - 1 for t > n . 


Adjoined Uni -Directional Half-Planes 
Figure 4 shows two uni -directional half-planes of different fiber 
and matrix properties which are assembled to form adjoined half-planes. 
Both planes may have an arbitrary number of broken fibers. Superscripts 
I and II are used to designate quantities corresponding to plane I and II 
respectively. The normalized spatial variables, n and 5 , in the longi- 
tudinal direction are related by 



The shear stresses, rjln) and at the interface are normalized 

a a 

with respect to material properties of plane I and II respectively, and 
are related to the actual shear stresses as follows: 
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Figure 4. 


(23) 


r. 


Tg (5) = 


1 


Epth^^ T^(y) 

*f'"h / 


II 


where, Ti(y) and x;^(y) are the actual shear stresses at the interface 

a a 

acting on plane I and II. Since the remote displacements and hence the 
remote strains are assumed equal, the remote stresses in planes I and II 
must be related by 
II 

(24) 


I 

a 


• I 


II 


- 1 - 1 1 

If T_(n) and t, (?) are known then the solution for each of the half- 

a a 

planes is the same as that developed in the previous section. Equation (.20) 
The Fourier constants and corresponding to broken fibers in plane I 
and II can then be obtained by solving the following sets of equations 
for known tJ(ti) and tJ^(?): 

Q a 


E / cos[(N*-j + m + ^)e]cos[(n + i) 0 ] 6 de 

m=l ^ 0 

pTT 1 oo~5tT 

+ -/ cos( 0 / 2 )cos[(n + -^) 0 ]/ e x^(t)dt d 0 = 1 , 

^ ^ a 


and E B|J^ / cos[(N *2 + m + -^)o]cos[(£ + i)0]6d0 , 
m=l 0 

p 7T 1 00 “(St TT 

+ -/ cos(0/2)cos[(£ + ^)0]; e T^(t)dtd0 = l, (25) 


for n = N* + 1, N* + 2,...,N^ and £ = + 1 , Ng + , 

where, M-j and M 2 are the number of broken fibers in planes I and II. The 
normalized displacements of a fiber in plane I and II are then given by 


Vj(n) 


7T 


7T 1 6n 

/ E B^ cos[(N? + m + i)0]cos[(n + i)0]e do 
0 m=l ^ ^ 

- I / cos[(n + i) 0 ] r D(6,n,t) 


xf(t)dt do 

a 


(26) 
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and v“(5) 
where 


^2 “ ^ 5 

^ z cos[(N2 + m + -^)0]cos[(n + -^)e]e de 

- 1/ cos . (p/ 2l cos [(n + 1)9] r D(6,^,s) ^”(s)ds de, 
■^0 ^0 ^ 


v;^(y). 



oo n 


and 



(y) 



II,, II 

oo n 


(5) . 


( 27 ) 


(28) 


When the above two half-planes are joined together the shear stress 

along the interface is unknown, but from equilibrium the shear stresses 
I T T 

Tl(y) and t (y) acting on each of these two half-planes must be equal 
and opposite. Further, as the shear stress is directly related to the 
distortion of the matrix from the shear-lag assumption, it follows that 
these stresses must be proportional to the difference in the displacement 
of the first fibers of plane I and II, These conditions result in the 
following two equations: 

x'(y) = - T“(y). and (29) 


tj(y) = (G„/h)’ 


v*(y) - 


v”(y) 


(30) 


where, (6jyj/h) is the equivalent shear stiffness of the interface. It is 
interesting to note that in a continuous elasticity solution one would 
match surface tractions and displacements at the interface while in the 
present discrete modeling solution the shear stresses are required to be 
equal and the shear-lag relation. Equation (30), takes the place of the 
displacement equality. 
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Substituting for the actual displacements in terms of normalized values 
using Equation (28), and recalling the relation between the actual and 
normalized shear stresses from Equation (23), Equations (29) and (30) re- 
duce to 




-I 

■'a 


(n) = vj(n) - 


(n)/R-| , and 

(G„/h)V(G„/h)^ 


(31) 

(32) 


where, 



1 






ApEph 


(33) 


Using Equations (26) and (27), V^(n) and vjj^(5) can be found and 
substituted into Equation (32) thus resulting in an integral equation 
for T?(n) in terms of the Fourier constants, and , and the nor- 
malized shear stresses, T?(n) and tJ^(c). The spatial variable c and 

a a 

the normalized shear stress t^^(c) may be eliminated from the above in- 

a 

tegral equation and also from the equations for the Fourier constants. 
Equation (25), resulting in the following set of governing equations; 
two series equations and one linear integral equation: 

Ml 

5: B^ - / cos[(Nt + m + i)e]cos[(n+l)e]6d6 
m=l 0 ^ ^ 

9 TT 1 °° ”6t _-r 

+ — / cos(e/2)cos[(n + -i-)0] / e TMt)dt de = 1 , (34) 

0 ^ 0 ^ 

^1 

T T ? ^ * 1 1 

ZB -f cos[(N« + m + -^)e]cos[(ji + :^)0]6d0 

m=l m 0 

— cos(0/2)cos[(£ + l)0] / e ^ xht)dt d0 = 1 , (35) 

0 ^0 ^ 


18 



for n = N* + and Jl = N 2 + 1 ,...,N 2 9 and 

M, 


T 9 ^ T * 1 "5ri 

i(n) = G-, - / E cos[(Nn +m + -^)0]cos(e/2)e de 
® ‘ ” 0 ro=l 


-I, 

T. 


71 ”2 II 


l,.n -^1'' 


■ ’^7'^it r 2: B cos[(N 2 + tn + -i-)e]e 
0 m=l 


de 


0 0 


2 cos^(e/2) 

6 


S,--i 

D(6,Ti,t) 


where, = (G„/h)’(h/G„)^ , G,g = fi, = 6/R^ . 


de Tg(t)dt 

(36) 


The Buffer Strip Laminate 

Figure 5 shows a finite width strip (buffer strip) between two half- 
planes of different or of the same fiber and matrix properties. One of 
the half-planes, plane I, and the finite strip, region II, may contain an 
arbitrary number of broken fibers. The normalized spatial variables n, C 
and c are related to each other in the same way as in Equation (22). The 
normalized shear stresses x^^(c), ^ 5 ^^(c) air's related 

to their corresponding actual shear stresses in the same manner as in 
Equation (23). Further, the remote strains in all three regions are 
assumed equal. 

The solutions for planes I and III can be obtained for known T^(n) 

a 

and tJ^^(c) from the half-plane solution and therefore we need to de- 
termine the solution for the finite strip, region II, only. This solu- 
tion is developed by considering the half-plane shown in Figure 3 with 
the following special condition. It may be assumed that the matrix be- 
tween the fiber Nw and (Nw + 1) splits all the way to infinity and that 
an external shear stress ^^^(y) is applied on the surface of the split 
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Figure 5. The Three Regions of the Buffer Strip Laminate 



as shown 1n Figure 6. The introduction of a split with or without the 
external applied shear stress results in two special equilibrium equa- 
tions for the two fibers on either side of the split [9, 10, 11] in 
addition to those valid for the free- edge fiber and the generic fiber, 
n, given by Equations (8) and (9). These equilibrium equations in the 
normalized form are as follows: 

^2yll 

_2- + = ^”(5) . for fiber 0 , (37) 

<l£ 

= 0 , for fiber n , (38) 

^ - ''nw = - ^b 

d^ ^ 

d^V^^ 

''nw+ 2 - C+l " Nw + 1 . (40) 

The left hand sides of Equation (39) and (40) can be reduced to the 
standard form by adding or subtracting a term * V^^^) on both sides. 

Making use of the transforms similar to Equation (10) and (11) and fol- 
lowing the procedure of the half-plane problem, the single equation 
specifying the transformed normalized displacement V^^(^,e) can be 
obtained as 

= ?“(5)cos(e/2) + [g(s) - ?P(5)]F^ , (41) 

dC 

where, 

9(5) = - v” , and = cos[(Nw + ^)e] - cos[(Nw + |-)e] . 

The solution to Equation (41) satisfying vanishing stresses and dis- 
placements at infinity and unit compression on the crack surface is then 
given by: 
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applied shear stress 



of the split 


Figure 6. Half- Plane with Matrix Split 
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^2 

z - / cos[(No + rn + i)0]cos[(n + i)0]d9 

I" ‘H’ ^ ^ t <- 

m=l 0 

« TT 1 °° -TT 

+ — / cos(0/2)cos[(n + ^)0] f e x“(t)dt d0 

+ “ cos[(n + i)0] / e /g(s) - ds d0 = 1 , (42) 

■^0 ^ 0 I ^ J 

for all broken fibers. 

If g(t) and are known, Equation (42) reduces to a set of 

linear algebraic equations in and can be solved directly. 

VJhen the above finite strip, region II, is introduced between the 

two half-planes, planes I and III, as shown in Figure 5, it results in 

two interfaces each similar to the one discussed in the problem of adjoined 

half-planes. The interfacial shear stresses x^(n), and 

a a D 

can be obtained in the same manner as in the previous solution, 
b 

that is, using the following relations: 

Tg(y) = - and 

Tg(y) = (G|v,/h)''^ between planes I and II, (43) 

and 

x“(y) =xIII(y). and 

T?^(y) =-(G»/h)^^ [v^^(y) - v^^^(y)] between planes II and 111,(44) 
b Nw 0 

where (Gj^^/h)^^ and (Gj^/h)^^ are the shear stiffnesses of the respective 
interfaces. By definition g(c) is given by 

- < • (45) 

Equations (43) - (45) along with the two equations for the Fourier con- 
stants from the stress boundary conditions on the broken fibers in region 
I and II, Equations (34) and (42), represent the complete solution for the 
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finite width buffer strip problem. In Equation (45) it must be noted 
that is the normalized displacement of fiber Nw of the finite strip, 
whereas is the normalized displacement of the first fiber of a new 

half-plane having applied shear stresses along the free edge and 

fiber/matrix properties of the finite strip. 

As in the previous problem, first the actual displacements are sub- 
stituted in terms of their corresponding normalized displacements, then 
the normalized displacements are substituted in terms of the normalized 
shear stresses and the Fourier constants. Further, and 

can be eliminated knowing the relationship between the actual shear 
stresses. Equations (43) and (44). Then the final set of governing equa- 
tions for the finite width buffer strip problem can be listed as follows: 
Ml 

i cos[(N? + m + i)e]cos[(n + i)e]6d0 

T m TT ^ I ^ 

m=l 0 

—fit 

+ -/ cos(e/2)cos[(n + i)e] / e x^(t)dt de = 1 , (46) 


2 IT 

1 I / cos[(N2 + m + ^)e]cos[(5, + i)e]6de 


m=l 


— /"cos(0/2)cos[(£ + ^)0] / e ‘ T^(t)dt d0 


-Sit 


-I, 


+ I / F^cos[(£+^)0] / e fg(s) - ^P(s)| ds d0 = 1 , (47) 

^0 0 Iw ^ J 

for n = + 1,...,N-| and + 1 , . . . ,N 2 > 

Ml 

^ ^il f cos[(N* + m + l)0]cos(e/2)e de 
m=l 0 

Mo X 

2 tT 9 * 1 "°1^ 

z G^-jR-j — / cos[(N 2 + m + ^)e]cos(0/2)e de 
m=l ^ 0 
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oo G TT 2 

_ f _li f COS [D(6,n,t) + G-J2 D(6^,n,t)]de ^^Ct)dt 

0^0 


^ D(5,n/R^,s)d0 [g(s) - T^^Cs)]ds, 

0 "^0 (48) 

g(?) = - i 7 cos[(N 2 + m + -^) 0 ]cos[(Nw + ^)e]e de 
m=1 ^ 0 


^ 1/ cos(b/ 2) cos[(Nw + l)e] D(6,t/R-. ,^)de ^^(t)dt 
r^2 TT 0 5 2 i a 

« T 7T cos[(Nw + i)e] « 

+ / i D(6,c,s)de g(s)ds 

1 /f cos^O/Z) , F^cos[(Nw4)6]' 


0 R. 


0^0 


D( 6 , 5 ,s)de ?J^(s)ds, 


(49) 


and 


^2 •• 6 5 


_ *^12^12 1 Z’ cos , (^ . cos[(Nw + 1)8] D(6,t/Ri ,g)de ?ht)dt 


n 2 TT 

0 R-j 0 


oo G.p TT cos[(Nw + i)e] 2 X , V 

+ / — / J ^ F D(6,C,s)d0 g(s)ds 


00 G . 1 TT 

- f f 


^cos!^0(6„,c,s) 

L Rj 

COs[(NW + -^)0] 9 

+ 5 — ^ F‘^D(5,^,s) 


d0 xj^(s)ds (50) 


where 


Rl = 


«2 = 
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I 


«1 = «/R-| , 

^2 ~ ^^^2 * 

®23 = • 

The solution reduces to two series equations and three linear integral 
equations . 

Uni-Directional Finite Width Strip with Broken Fibers 
In this section a solution is obtained for a uni -directional finite 
width strip containing an arbitrary number of broken fibers as a special 
case of the solutions already developed. With reference to Figure 5 it 
can be seen that if the interfacial shear stresses (y) and (y) 
acting on the region II are zero then it results in a finite width strip 
with broken fibers subjected to a uniform remote axial stress as shown 
in Figure 7. Therefore the solution to the finite width strip problem 
can be obtained as a special case of the solution developed for region 
II of the buffer strip laminate by setting and tJ^(s) to zero in 

Equations (42). With reference to Figure 7 these equations are given as 

i cos[(N* + m + ^)e]cos[(n + ^)e]de 

m=l 0 

o o 1 °° -5t 

+ - / F^ cos[(n + i)0] / e g(t)dt de = 1 , (51) 

^0 0 

for all broken fibers, and 

M p TT 1 

g(^) = - z ^m “ cos[(N + m + ^)0]cos[(Nw + ^)e] e d 0 
m = 1 ^0 

■] 

CO « IT cos[(Nw + -h-)0] 9 

+ / ^/ ^ F'^ D(6,e,t)d0 g(t)dt (52) 
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Figure 7. Finite Width Uni-Directional Strip 
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where 


= cos[(Nw + ^)e] - cos[(Nw + f)e] » and 

g(5) = . 

Since by definition is the displacement of the first fiber of 

a new half-plane having fiber/matrix properties of the finite width 
strip and no applied shear stress along the free edge, it must be equal 
to zero. Therefore g(^) is given by 

g(c) = - . (53) 

Therefore the solution of the finite width strip reduces to one 
series equation coupled with one linear integral equation. The location 
and the number of broken fibers are arbitrary except that the fiber 
breaks must be along the x-axis. The solution obtained holds for a 
central notch, an edge notch, an off-center notch, or for multiple 
notches along the x-axis. 
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CHAPTER III 


SOLUTION TECHNIQUE 


In all the problems dealt with so far, the solution reduces to one 
or two series equations coupled with one or more linear Fredholm in- 
tegral equations of the second kind. Since there is no exact closed 
form solution available to solve such a system of equations, a numerical 
procedure is developed and presented in this section. The technique 
makes use of a method given by Riez [12], to solve a linear Fredholm 
integral equation of the second kind defined within a semi -infinite 
interval of integration. The solution is based on the fact that a given 
integral over a semi-infinite interval may be approximated by the 
Gauss-Laguerre quadrature rule as 
00 k -x^- 

/ f(x)dx = Z w. e f(x. ) (54) 

o i=l ^ 

where x^- is the i"^^ zero of the Laguerre polynomial, L|^(x), and w^- is 
the corresponding weight function given by 

w. = 2 . (55) 

' [(k+1) Lk+i(xi)]‘^ 

The Laguerre polynomial L|^(x) is given by 


. (X) = . 


(56) 


dx’ 


Since the form of the equations for each of the solutions is the 
same, the development and application of the numerical procedure can be 
demonstrated, without loss of generality, by taking the equations corre- 
sponding to one of the solutions. Consider then the solution corre- 
sponding to the problem of adjoined half-planes given by Equations 
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> T o °° 1 1 

j; Bz — / cos[(N? + m + 4)e]cos[(n + i)e]sde 
1=1 0 ' 

? TT r- ^ “5t _T 

+ — / cos(e/2)cos[(n + :^}0] / e TMt)dt de 

TT c ^ a 


E / cos[(N 2 + m + ^)e]cos[(£ + -^)e] 6 de 

m=l 0 

- — ^ / cos(e/ 2 }cos[(£ + i)e] / e ^ T^(t)dt de = 1 , ( 58 ) 

for n = N*+ 1 ,...,N-| and £ = + 1 , . . . ,N 2 > and 

TT ^1 

T^(n) = ^ cos[(N* + m + l)e]cos(e/ 2 )e de 

1^2 -5 ^ 

- i ; E C 0 s[(N 2 + m+ ^)e]GOs(0/2)e ^ de 

^ 0 m=l 


■ Si D(6,n,t) + ^ D(6pn,t) 

0 0 1 

X de ?i(t)dt . (59) 

a 

Integrals over the spatial variable, t, in the above equations, 
defined over the semi -infinite interval, can now be replaced by the above 
series representation to yield 


^1 f TT 1 

E^ < I ; cos[(N* + m + ^)e]cos[(n + ^)e] 6 de 

k C n -K 1 - 6 t-| I 

+ E | — / cos(e/ 2 )cos[(n + i)e]e de w-e 

j=l 0 

Sr , 1 

E (“/ cos[(No + m + -^)e]cos[(£ + 4 -)e] 6 de > B", 

1=1 0 J 
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k 

- z 
j=l 


Gip O IT 1 

-S-f cos(e/2)cos[(5, + y)e]e de w.e 


( 61 ) 


for n = N* + 1 , . . . ,N-| , and A = + 1 , . . . »N 2 , 


and 




+ z K(n. , tjw.e 
ij 1 J J 

^1 f p TT*^ * 1 “*^^1 1 I 

- z ( G--| ^ / cos[(N-| + m + -^)e]cos(e/2)e 
m=l 1^0 J 


Mp ^ 


+ z 
m=1 


-6^n^ 


R,G^^ ^ cos[(N 2 + m + ^)0]cos(e/2)e ' ‘ de 

0 

X = 0, i = 1 k 

where, S.. = 1 for i = j 


(62) 


ij 


=0 for i ^ 3 > 

and K(nptj)= [D(6.n^.tj) + -^ D(6-| ,n,- ,tj)]de . 

Therefore, the solution of the integral equation coupled with a set of 
series equations reduces to solving a system of linear algebraic equations 
in and explicit values of T^(n) at the quadrature points. The 

above system of equations is solved by the method of Gauss-elimination 
with partial pivoting. 

It must be pointed out that the terms in the integral equation of 
the form have cusp at n = t and the quadrature representation in 

the standard form. Equation (54), results in large inaccuracies especially 
for large values of n. This has been overcome by modifying the above 
quadrature rule by deleting selected terms in the series and replacing 
them by closed form integration in the neighborhood of the cusp. With 
reference to Figure 8, the integral over the semi -infinite interval is 
then represented by the modified equation 
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/” KCn^.t) ^g(t)dt = z K(n^,t-) ^ w. 

o j=l 

[ t • t- 


-I 




1 


/ K(n..t) ^J(t)dt+ / K(n.,t) ? 3 (t)dt 

1 a 1 a 




(63) 


where. 


= (t._i - t.^ 2^/2 . and 

^2 ^ ^^ i+2 " ■ 


Since x^(t) is yet unknown, the two integrals on the right hand side 
a 


of Equation (63) can not be evaluated in a closed form. However, if the 

function T^(t) is taken to be the average of the terminal values within 
a 

each of the intervals, the above two integrals can be expressed as 




/ K(n.,t) *J(t)dt = — 5 2 ^ K(ti-,t)dt 


and 


/ K(rn.t) xht)dt = 2"-^ / K(nT,t)dt . (64) 

t. t. 

Substituting Equation (64) into Equation (63) and rewriting 
results in 

-I 


-I,, 


/ K(n^.t) Tg{t)dt = z* K(n^.t) ^a(tj)e '' Wj 
o j-1 

1 I 

+ ^ / K(n^.t)dt Tg(t^.i) 
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t.^^+A2 


f K{n^.t)dt+ / K(ri^-,t)dt 


S‘-T*l 




*1+1 "^^2 




+ ^ ' } ‘K(ni,t)dt . 


( 65 ) 


*i 


where E exclu(ies terms corresponding to j = i - 1 , i and i + 1 . 

The quadrature-rule in the modified form. Equation (65), is used to 
represent those integrals in which the integrand has terms with cusp. 

The two series equations remain unaltered and the substitution of Equa- 
tion (65) into Equation (59) results in 


s*K(ni,tj)Wj e + a._-j ) + U+a- )yt.) + a-+^ ^a^^i+1^ 

0 • 

f 2 * 1 "“^^i 1 I 

- E cos[(N-j + m + ^)0]cos(0/2)e f 

m=l 1^ ^ 0 J 

M. r 

^ ^ * 1 
+ E cos[(N 2 + m + ^)0]cos(0/2)e 

m=l I ^0 

^ IT. 


d0 


X = 0 , 1 = 1 


where, 


*i 


i-1 =K ^ K(ni,t)dt . 


*i-T*l 


t ^+^+&2 


I. =4- / K(m,t)dt + -i- / K(n.,t)dt . and 

T ^ 4. A ’ ^ + 

*1 


*i+r^^2 


‘i+1 2 


= i / K(n,-,t)dt 


*i 


This again reduces to solving set of linear algebraic equations as before. 
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CHAPTER IV 


RESULTS AND CONCLUSIONS 

First, typical results are given for a Graphite/epoxy laminate con- 
taining a buffer strip as shown in Figure 9. The crack growth behavior 
of the lamina is studied by varying the buffer strip material, the width 
and the thickness of the strip and the initial crack length. The three 
materials considered for the buffer strip are Nylon, Kevlar and S-glass, 
all in an epoxy matrix. The matrix and fibers are assumed to be linearly 
elastic. The failure criterion is that a fiber fails upon reaching its 
ultimate failure stress as determined from an unnotched coupon test. 

Figure 10 presents results corresponding to initial crack growth in 
plane I, crack arrest at the interface, crack growth in the buffer strip 
and subsequent lamina failure. In these results all fibers are of the 
same cross-sectional area and in all cases the buffer strip is ten 
fibers wide. Since all the buffer strip candidate materials are of 
lower modulus than that of the parent laminate, which in this case is 
Graphite/epoxy, the stress concentration factor at the near end of the 
notch (nearer to the interface) is always higher than that at the far 
end of the notch. The solid line in Figure 10 represents the remote 
stress required to initiate crack extension, (fail the first unbroken 
fiber in front of the notch, fiber A in Figure 9). The remote stress 
required to fail the lamina catastrophically, (fail the first fiber in 
plane III, fiber B in Figure 9) is given by the broken line in Figure 10. 
Both these stresses are functions of the initial crack length and decrease 
with increasing length. The crack growth takes place by breaking 
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consecutive fibers from the crack tip to the interface. Then, depending on 
the stress level required to run the crack to the interface and depending 
on the buffer strip material, the crack may arrest. It is very interesting 
to note that all three buffer strip materials require an increasing stress 
to continue the crack growth in the buffer strip, although Kevlar will 
only arrest a crack if it initiates under fairly low load, i.e. initially 
close to the interface. 

Total lamina failure will occur when either the buffer strip is fully 
broken or when the first fiber in plane III, fiber B, fails. In both 
cases continued crack growth is unstable once fiber B fails. For the 
particular lamina of Figure 9, all fibers in the Kevlar buffer strip fail 
before fiber B attains its failure stress, whereas for S-glass and Nylon, 
fiber B fails when there are still some fibers left unbroken, i.e., the 
crack jumps the buffer strip. The ultimate lamina failure stresses are 

cf^ = 0.272 for Kevlar 

= 0.395 for S-glass, and 

= 0.444 for Nylon, 

where is the ultimate fiber stress as determined from an unnotched 

Graphite/epoxy laminate. The material properties used for these results 
are given in Table 1. The results of Figure 10 indicate Nylon to be the 
best of the three materials but this is only true if the matrix can 
support the very large failure strains (about 20%) of Nylon. In a typical 
angle-ply laminate with Nylon buffer strips and with continuous +45 
graphite plies, high strain levels certainly cannot be reached before 
failing the angle-plies and continuing the crack. This behavior was 
observed by Verrette and Labor [4]. The extension of this study to 
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account for angle plies is surely necessary to represent accurately such 
low modulus fibers. 



TABLE 1. 

FIBER 

PROPERTIES 


f i ber 

mass dens 
Gms/m^ 

ity 

modulus 

MPa 

ult. stress 
MPa 

Nylon 

84 


2000 

233 

Kevlar 

106 


111400 

2020 

S-glass 

199 


101000 

2800 

Graphite 

155 


300000 

2800 


In Figures 11-13 the effect of buffer strip width on crack growth 
through the strip is indicated. The ultimate failure stress of the 
lamina as a function of buffer strip width is plotted in Figure 14. 

From Figure 14 it is seen that for Nylon the ultimate failure stress 
continues to increase with an increase in width, whereas for S-glass and 
Kevlar about ten fiber widths is optimum. 

Table 2 gives a comparison of the ultimate lamina failure stress 
for a buffer strip of ten fibers, first for equal area and then for 
equal weights of the fibers. These results are normalized with respect 
to S-glass. A Kevlar region of two plies is approximately equivalent to 
one ply of S-glass in weight and it is seen that S-glass is still superior 
to Kevlar. This agrees with results obtained by Poe and Kennedy [5]. 

Next, results are presented for the finite width strip with broken 
fibers forming a central notch and subjected to a uniform remote axial 
stress as shown in Figure 7. The stress concentration factors for various 
notch lengths (number of broken fibers, M) are obtained corresponding to 
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Figure 11. Effect of Buffer Strip Width 
on Crack Growth ( Kevlar) 



Figure 12. Effect of Buffer Strip Width 
on Crack Growth (S-Glass) 




Figure 13. Effect of Buffer Strip Width 
on Crack Growth (Nylon) 


Figure 14. Ultimate Lamina Failure Stress 
Vs. Buffer Strip Width 



TABLE 2. COMPARISON OF LAMINA FAILURE STRESS 


fiber 

O' /cr..n. 
equal area 

^w^*^Ul t 
equal weight 

Nylon 

0.444 

0.472 

Kevlar 

0.272 

0.318 

S-glass 

0.395 

0.395 


two widths, Nw = 24 and 48. Following Hedgepeth [6] the stress concentra- 
tion factors for a uni -directional infinite region are obtained using the 
relation 


where 


°N+1 _ 4-6-8 (2H+2) 

'’N+I ■ 3-5-7 (2M+1) 


(67) 


N+1 = index of the first unbroken fiber at the notch tip, 
M = total number of broken fibers in the notch, 

^N+1 " stress concentration in fiber N+1, 

Oj^^l = axial stress in fiber N+1, and 
= applied remote stress. 


By comparing the above two stress concentration factors corresponding 
to a given notch length the finite width correction factor for a uni- 
directional strip is obtained and is given in Table 3. Also given in 
Table 3 are the corresponding finite width correction factors for an 
isotropic strip obtained from the following relation [13] 

Yj = 1.0 + 0.1282 (2a/W) - 0.2881 (2a/W)^ + 1 .5254 (2a/W)^, (68) 

where, 

'^2 = finite width correction factor, 
a = half-crack length, and 
W = width of the strip. 
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From Table 3 it is seen that the finite width correction factors 
for the isotropic strip and the uni -directional strip are not signifi- 
cantly different for aspect ratios (notch length to strip width ratios) 
less than 0.5. For higher values of aspect ratios there is a considerable 
difference between the two cases. An aspect ratio of 0.9167 corresponds 
to 22 broken fibers in the 24 fiber strip and all the applied load is 
carried by the two remaining intact fibers giving a stress concentration 
factor of 11.0, (22x1.0/2). Similarly, the other limiting case corre- 
sponding to the aspect ratio of 0.9583 in which 46 fibers are broken out 
of 48, leaving the two end fibers to carry all the applied load giving a 
stress concentration of 23.0 (46/2). These two limiting cases are pre- 
dicted accurately by the above solution indicating the validity of the 
model . 


TABLE 3. FINITE WIDTH CORRECTION FACTORS 






Fini te 

width correction factors 





2a/W 

uni -directional strip 

Isotropic 



a — ^ 
[= - •■=» 



Nw = 24 

Nw = 48 

1 

1 

0.1667 

1.01334 

1.01125 

1.02042 

1 



1 

0.25 

1 .02890 

1.02514 

1.03788 

■ 

^ U 

1 

0.3333 

1.05207 

1 .04581 

1.06722 

1 



1 

1 





0.5 

1.13040 

1.11598 

1.18275 

2a = M X 

d 

0.75 

1.44760 

1 .39280 

1.57762 

W = Nw X 

d 

0.9167 

2.57412 

3.7 

2.05049 

c 

l=fiber 

spacing 

0.9583 


3.77560 

2.20070 
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From the above results and discussions the following conclusions 
can be made. The method predicts the fracture behavior of a hybrid 
laminate in terms of material properties, geometry and initial crack 
length. The results agree well with those obtained experimentally for 
buffer strips of high modulus fibers such as S-glass and Kevlar, where 
the stiffness of the angle-plies is very small compared to that of zero 
degree ply. For low modulus (high failure strain) buffer strip materials 
such as Nylon, in which the stiffness of continuous angle-plies is com- 
parable to that of the axial buffer strip fibers, the model is inadequate 
due to the basic assumption made in idealizing the laminate as a uni- 
directional composite. The method predicts the best buffer strip material 
to be one with a low modulus and as large an ultimate strength as 
possible, e.g. , S-glass is superior to Kevlar. 

As stated earlier, the main aim of this work was to understand the 
basic mechanism of crack growth and arrest in hybrid laminates, keeping 
the model as simple as possible. This goal has been achieved. However, 
in order to represent an actual buffer strip panel more realistically, 
the model certainly needs and has the potential for, further modifications. 
An immediate extension is to model the panel as a periodic uni-directional 
hybrid laminate simulating the regular placement of the buffer strips. 

Next, the role of angle plies must be accounted for by either adding a cover 
sheet over the uni -di recti onal laminate or by some other means. Finally, 
the effect of matrix damage in the form of longitudinal yielding and 
splitting at the crack tip and at the interfaces must be incorporated. 
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unidirectional finite width strip containing an arbitrary number of broken fibers. 
The stress concentration factors for the finite width strip are compared with 
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